Creation of a dipolar superfluid in optical lattices 
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We show that by loading a Bose-Einstein condensate (BEC) of two different atomic species into an 
optical lattice, it is possible to achieve a Mott-insulator phase with exactly one atom of each species 
per lattice site. A subsequent photo-association leads to the formation of one heteronuclear molecule 
with a large electric dipole moment, at each lattice site. The melting of such dipolar Mott-insulator 
creates a dipolar superfluid, and eventually a dipolar molecular BEC. 



The physics of strongly correlated systems in ultra- 
cold trapped quantum gases has attracted a growing in- 
terest recently. The experimental progress on the trap- 
ping and cooling techniques, and the control of the inter- 
atomic potentials by means of Feshbach resonances allow 
to analyze situations beyond the validity of the mean- 
field approach. In particular, the possibility to achieve 
ID bosonic gases opens the pathways toward a gas of im- 
penetrable bosons [Q. Also, a rapidly rotating trapped 
bosonic gas should exhibit effects similar to the frac- 
tional quantum-Hall effect ||. The field of strongly 
correlated atomic systems concerns also ultra-cold Bose 
gases with large scattering lengths and the Mott- 

insulator (MI) to superfluid (SF) phase transition || in 
cold bosonic gases in optical lattices ||F?j] . The latter has 
been recently observed in the remarkable experiment Q . 

The physics of ultra-cold molecular gases is also a very 
active research area. Cooling and trapping of molecules 
can be achieved by means of buffer-gas techniques || , by 
employing the dipolar moments of polar molecules pO[ , 
and by means of photo-association of already existing 
ultra-cold atomic gases 11 . Jaksch et al. have re- 



cently proposed that a bosonic atomic gas, which is in a 
MI phase with exactly two atoms per site in an optical 
lattice, can be used to obtain via photo-association a su- 
perfluid gas of homonuclear molecules, and eventually a 
molecular BEC. 

The influence of dipole-dipole forces on the properties 
of ultra-cold gases has also drawn a considerable atten- 
tion. It has been shown that these forces significantly 
modify the ground state and excitations of trapped dipo- 
lar BEC's [p~3| — P~S|] - In addition, since dipole-dipole inter- 
actions can be quite strong relative to the short-range 
(contact) interactions, dipolar particles are considered to 
be promising candidates for the implementation of fast 
and robust quantum-computing schemes |19 20[ ). Several 
sources of cold dipolar bosons have been proposed, in- 
cluding atoms with large magnetic moments Q], and 
externally induced electric dipoles E3Jla,tt9]. Perhaps 



the most promising perspective in this sense is provided 
by ultra-cold gases of polar molecules with large electric 
dipole moments, which could eventually dominate the 
physics of the system ]l5],[l6| . 

This Letter is devoted to the analysis of the genera- 
tion of ultra-cold polar molecules in an optical lattice. 
We consider the loading of the lattice by a superfluid of 
two different atomic species, and analyze the transition 
into a MI phase with exactly one atom of each species per 
site. This transition is followed by the formation of polar 
dimers by photo-association, or adiabatic passage. Fi- 
nally, the quantum melting of the molecular Mott phase 
into a polar molecular superfluid takes place. 

We consider a miscible gas of two trapped atomic 
bosonic species at very low temperature T. Although 
our calculations could in principle be employed for arbi- 
trary atoms, we have assumed in the calculations below 
the particular mixture of 41 K and 87 Rb, whose simulta- 
neous Bose-Einstein condensation has been recently re- 
alized [^lj . We additionally consider that the binary gas 
is confined in a 2D optical lattice, similarly to the 3D 
one-component case recently observed in Ref. The 
transverse dimension is considered as harmonically con- 
fined. The physics of the system is governed by the two- 
species generalization of the Bose-Hubbard (BH) model, 
described by the Hamiltonian: 



H = 2_j I J a a i a j + J bb\bj] + Uab ^ U ai n bt 
<i,j> i 

+ ^ ^[Uoariaiinai - 1) + U Qb n bi (n M - 1)], (1) 



where a^, bi are the annihilation operators of 41 K and 
87 Rb atoms at the lattice site i, which occupy a state 
described by the Wannier functions w a (r — r^), w b (r — 
rj) of the lowest energy band, localized on this site. 
This implies that energies involved in the system are 
small compared to the excitation energies to the sec- 
ond band. We denote the position of the local min- 
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imum of the optical potential as r,-, and the number 
operator of K (Rb) atoms at the site i as n a i = a\ai 
(riw = b\bi). In Eq. (|l|) only the nearest-neighbor tun- 
neling is considered, which is described by the parameter 

Ja,b = /< 6 (r - rOhfeV 3 + Vz(r)]«; a , 6 (r - r,)d 3 r , 
where j and i are the indices of neighboring sites, and 
^( r ) = J2^=x y cos 2 (fc^^) is the optical lattice po- 
tential (with wavevector k) which we consider to be 
the same for both species. The optical potential for 
K and Rb is almost the same for the detunings corre- 
sponding to the wavelength 1064nm of a Nd:Yag laser 
p2f . The same species interact via a short range pseudo- 
potential, which is described by the coefficients Uo a , Uob'- 
U 0a ,b = 4-7r^ 2 a K ,Rb J |wa,6(r-ri)| 4 /m K ,Rb d 3 r, where ok, 
wik (ctRbj ^Rb) are the scattering length and atomic 
mass of potassium (rubidium). The interactions be- 
tween different species are characterized by the coefficient 
U ab = 2irh 2 a KRh f \w a {r - r,-)| 2 |tu 6 (r - r, i )| 2 /m ab d 3 r, 
where olkko is the K — Rb scattering length, and m Q f, 
is the reduced mass for K and Rb atoms. The scatter- 
ing length aKRb has been estimated recently to be of the 
order of 8.6nm j2l), that is the repulsive interactions be- 
tween K and Rb are stronger than those between K-K 
and Rb-Rb, which are characterized by ok = 3.17nm 
and ORb = 5.24nm, respectively. As a result, it is ener- 
getically favorable for K and Rb particles to remain im- 
miscible, i.e. to stay at different lattice sites, rather than 
in the same ones. This situation can be however modi- 
fied by means of Feshbach resonances, which can strongly 
modify the value of OKRb- We have thus assumed that 
the value of aKRb can be reduced to 0.5y'aKaRb- In that 
case the two components acquire a miscible phase. 

In the 2D calculations presented below, we have con- 
sidered square lattices up to 80 x 80 sites, with periodic 
boundary conditions. In the first stage, we demonstrate 
the possibility of performing a transition from a state 
with two miscible superfluids to a Mott insulator phase 
with exactly one Rb and one K atom per site. To this 
aim we first obtain the ground state of the system using 
a variational Gutzwiller ansatz |jj . 

We approximate the ground-state wave function by 
\^mf) = Hi \4>i)i where the product is over all lattice 
sites. The functions \<pi) for each site are expressed in 
the Fock basis, {(pi) = J2n,m=o fn;m\n, m)i, where n, m 
indicate the occupation number for the K and Rb atoms, 
respectively. For numerical reasons we have assumed 
that maximally 3 atoms of each species can be present at 
any lattice site, and checked the validity of this assump- 
tion self-consistently. The coefficients {/n*L} are found 
by minimization of (^mf\H - YlilVaUai + jU&n&j] \^mf), 
where the chemical potentials \x a , fJ-b were chosen so that 
the mean numbers of K, Rb atoms were equal to the 
number of lattice sites. 

The ground state for weak lattice potentials corre- 
sponds to the coexistence of superfluids for both atomic 



species, and it was used as the initial condition for subse- 
quent time-dependent calculations. This state is charac- 
terized by a non-vanishing value of the superfluid order 
parameters \{a,i)\, \(bi)\ (Fig. |l|). To drive the system 
from superfluid to MI phase we slowly increase the lattice 
potential. This process can be described using a dynam- 
ical Gutzwiller method, as discussed in Ref . Q . To this 
aim one allows the coefficients {fn}m} to be time depen- 
dent and to follow the dynamics derived from the time 
dependent variational principle. The resulting equation 
becomes: 
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When the lattice potential increases, the system fol- 
lows the changes quasi-adiabatically, and the two species 
enter sequentially the Mott phase with one atom per site. 
As shown in Fig. |l| this occurs first for the heavier Rb 
atoms at t as 0.75s, and then at t ss 1.5s for K atoms. The 
ratio —Uo a /J a (—Uob/Jb) at the transition point equals 
approximately 23.5 (26), which interestingly compares 
very well with the expected value, 23.2, for a 2D single- 
component gas. Therefore the transition point is rather 
unaffected by the interactions between different species. 

At the end of the process each site contains one Rb 
and one K atom, with a relative atom number fluctu- 
ations less than 3%. This, however, requires a suffi- 
ciently slow modification of the lattice potential, typically 
within a time scale of 2 seconds for a lattice wavelength 
A = 1064nm. Note that the time scale is provided by 
the inverse of the recoil energy of the lattice, and there- 
fore it scales as A 2 . We note also at this point that a 
transition 3 times faster leads to about two times larger 
number fluctuation in the final state. We want to stress 
that during the whole time evolution the mean number of 
particles was constant, and equal to the number of lattice 
sites, which gives us the most complete analogy, within 
the Gutzwiller approach, to the experimental situation. 

Having obtained a Mott phase with one Rb and one 
K atom per site, a heteronuclear molecule can be gener- 
ated at each lattice site. One possible method would be 
two-color Raman photoassociation to make ground state 
1 E + or 3 E + dimcrs. Since heteronuclear species do not 
have gerade/ungerade symmetry, levels of either symme- 
try could be formed via intermediate excited states of 
fl = + or 1 symmetry, where Q labels projection of 
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electronic plus spin angular momentum on the molecular 
axis. Since the ft states excited at long range typically 
have mixed singlet and triplet character, photoassocia- 
tion to ground 1 S + levels is even possible for collision of 
two doubly spin-polarized atoms, which collide via the 
ground 3 £ + potential. The 1/R 6 dependence of the ex- 
cited molecular potential for the case of heteronuclear 
molecules, in contrast to the 1/R 3 dependence for the 
homonuclear case, leads to excitation at relatively short 
internuclear separation R. This difference may lead to re- 
duced Franck-Condon factors for the first exciation step 
but enhanced factors for the second de-excitation step 
. In any case, these can be worked out once the poten- 
tial energy curves of the KRb dimer are better known. A 
second possible method of making ground state molecules 
is available for heteronuclear dimers, namely, use of a 
pulse of microwave radiation to directly associate two col- 
liding atoms in a one-color transition from a scattering 
state to a bound dimer state. However, the dipole matrix 
elements within the 3 E + state are small, so this will need 
fairly high microwave power. Either method will result 
in an excited vibrational level not far below the disso- 
ciation limit. It should be possible to use a succession 
of Raman pulses to transfer the population from such a 
level to the ground v = vibrational level of the 1 E + 
or S + state, in a manner similar to that described by 
Rcf. pJ; microwave/infrared pulses may also be feasible 
for this heteronuclear species. The 3 £ + v = level has 
a small dipole moment of 0.0046 atomic units (1 atomic 
unit is 8.478 x 10 -30 Cm), whereas it is much larger, 0.30 
atomic units, for the 1 S + v = level In the latter 
case, special care should be taken to prevent the heating 
due to black-body radiation, which could be avoided by 
an appropriate shielding. 

After the creation of the K-Rb dimers, the system be- 
comes single-component. However, since the molecules 
have a large dipole moment and are considered to be ori- 
ented by an external electric field (« lOOV/cm), we have 
to include the dipole-dipole interactions in our calcula- 
tions. The corresponding Hamiltonian reads 0: 

<i,j> i 

<i,j> <<i,j» 

where hi is the annihilation operator of a heteronuclear 
molecule at the lattice site i. The number operator rij, 
the Wannier basis w(r — r,) and the tunneling coefficient 
J for molecules are defined similarly as in Eq. ([!]) . How- 
ever, the interaction part significantly differs from that 
of Eq. (0). The coefficients U a are defined as follows: 
U a = J \w(r-ri)\ 2 V mt (r—r')\w(r' — r 3 )\ 2 d 3 r dV, where 
a = |rj — rj||k|/47r is the dimensionless intersite distance. 
In particular, Uq determines the on-site interactions, U ai 
the nearest-neighbor interactions, U a2 the interactions 



between the next-nearest neighbors, etc. Consequently, 
the respective summations in Eq. (^) must be carried 
out over appropriate pairs of sites which are marked by 
< . . . > for the nearest neighbors, << ... >> for the 
next-nearest neighbors, etc. In the 2D calculations pre- 
sented below, we have taken into account interactions 
with up to 4 neighbors (<ji = 1, 02 = \/2, 03 = 2, 
(74 = v5), since the effects of interactions of a longer 
range are negligible. In the case of polarized dipoles the 
interaction potential is 

, 9 l-3cos 2 6> Anh 2 a c . . , 

V ^ = d k - r1 3 +-flr*fr- f, )» (4) 

where the first part of Eq. (|J) provides the dipole-dipole 
interaction characterized by the dipole moment d and 
the angle 8 between the dipole direction and the vector 
r — r'. The second part of Eq. (Q) represents the short- 
range interactions given by the s-wave scattering length 
a and the molecular mass M. We use the dipole moment 
as 0.3 a.u. |2j|, and assume the s-wave scattering length 
of K-Rb molecule as a = y / OKORb [ p5[ . 
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FIG. 1. Dynamical transition from the SF phase of 41 K 
(solid line) and 87 Rb (dashed line) atoms to the MI phase. 
The upper plot shows the value of the order parameters |(<ii}|, 
\{bi)\ (constant for all lattice sites) for both species, while 
the lower one depicts the variance Var(n) = U in 2 ) — In) 2 
of the on-site occupation. The lattice potential was changed 
as Vo(t) = Vsf + (Vmi - V SF )(t/t f exp(l - (i/t ) 3 ), where 
Vsf — 4, Vmi = 19 in units of recoil energy of 41 K atoms, 
to — 1.8s for the case of lattice wavelength 1064nm. 

We have checked that the MI phase of these heteronu- 
clear molecules is indeed the ground state of the system 
for a mean lattice filling factor equal to 1. We employ 
this state as the initial condition for our dynamical calcu- 
lations of the melting of the Mott phase towards a molec- 
ular superfluid. It is, however, interesting to note, that 
by itself the MI phase of dipolar particles offers promis- 
ing perspectives as a quantum computation device p0| . 
The melting dynamics is analyzed by employing again a 
dynamical Gutzwillcr ansatz, but with the Hamiltonian 
(Q). During the dynamics we consider a first stage in 
which we reduce the lattice trapping potential up to 4 in 
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recoil energy units of K — Rb molecules. This first stage 
allows us to reduce the value of — Uq/J, but the system 
remains still in the MI state mostly due to the large con- 
tribution of the repulsive dipolar on-site interactions to 
the Uo coefficient. This process can be performed in a 
relatively fast time scale compared to the total duration 
of the melting process, since the MI gap guarantees the 
adiabaticity at this stage. 
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FIG. 2. Melting of 41 K - 87 Rb dimers initially in the MI 
phase. The plot shows the time evolution of the molecular 
superfluid order parameter (solid line), which is the 

same for all lattice sites. The dashed line refers to static 
calculations of the ground state of the dipolar molecules 
placed in the lattice. Note a very sharp transition from 
the MI phase to the SF one, that occurs at t ~ 0.2s. 
The frequency of the transverse confinement was changed as 
uj z (t) = lomi + (u>sf - u>Mi)t/to, where lomi/2-k = 8kHz, 
u)sf/2it = 6kHz, and to = 0.65s. 



Later on, we have observed that the dynamics can 
be better controlled by modifying the transverse con- 
finement. As shown in Ref. J7||l5|], the properties of 
the dipolar gas strongly depend on the aspect ratio of 
the on-site wavefunction. In particular, by reducing 
the transverse confinement, the on-site wavefunction be- 
comes more elongated in transverse direction. As a con- 
sequence, Uq decreases M, and the system enters into the 
superfluid phase. Fig. Q shows this process. Due to the 
non-completely adiabatic evolution the gas does not ex- 
actly follow the expected stationary result, although very 
clearly a superfluid molecular phase is accomplished. The 
ratio — Uq/J 23.2 at which the transition from MI to 
superfluid equals surprisingly to the value expected for 
a 2D single-component gas without a permanent dipole 
moment. As shown in Fig. ^| this dynamical transition 
into the superfluid phase occurs for realistic parameters 
in a time scale of Is. This transition can be realized, 
50 times faster, although the dynamics of the superfluid 
order parameter is in such case not as smooth as in Fig. ^. 

In this Letter we have analyzed the formation of a 
superfluid of polar molecules. We have considered the 
loading of a BEC of 87 Rb and 41 K atomic species into 
an optical lattice. The modification of the lattice poten- 
tial produces the dynamical transition into a MI phase 



with only one atom of each species per site. The atoms 
can then be photo-associated on-site. Once the polar 
molecules are created, a molecular MI is formed. By re- 
ducing the lattice potential and modifying the transverse 
confinement, such MI evolves into a molecular superfluid 
within a feasible time scale. Eventually, a molecular BEC 
is formed. The mechanism studied in this Letter, pro- 
vides a perspective towards a dipolar BEC. 
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